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ABSTRACT 

We have mapped the [Ci] 3 Pi- 3 Po emission at 492 GHz toward the supernova 
remnant Cas A. We detect [C i] emission from the periphery of the diffuse Photon 
Dominated Region (PDR) covering the disk of Cas A, as traced by the carbon recom- 
bination lines, as well as from the denser PDRs associated with the molecular clouds 
towards the south-east. [Ci] emission is detected from both the Perseus and Orion 
arm molecular clouds, with the —47 km s _1 Perseus arm feature being strong enough 
to be detected at all positions. We estimate the C/CO relative abundance to be 0.2 
at the position of the identified CO clouds and > 1 for most of the cloud. Here we 
show that the distribution of [C i] emitting regions compared to the C + region and 
molecular cloud is consistent with a scenario involving PDRs. Using physical models 
for PDRs we constrain the physical properties of the [C i] line-forming regions. We 
estimate the densities of the [Ci] emitting regions to be between 10 2 and 10 3 cm~ 3 . 
Based on rather high volume filling factors (~ 50%) we conclude that [C i] emission 
mainly arises from diffuse neutral gas in the Perseus arm. 

Key words: ISM:general - ISM:lines and bands - ISM:molccules - ISMxlouds - 
radio linesdSM - radio lines:general - Galaxy:general. 



medium (diffuse PDRs). PDRs are defined as neutral regions 
where th e chemistry and heating are regulated by the FUV 
photons llHollenbach et al~I IT999). Current understanding 
of PDRs suggest chemical stratification in which with in- 
creasing depth from the surface of the PDR, the dominant 
carbon-bearing species changes from C + through C° to CO. 
Observations of [Cn], [Ci] and CO lines tracing the differ- 
ent layers of PDRs is thus a useful tool to constrain various 
conditions of the PDRs. 

Ionized carbon in dense star forming regions is mainly 
traced using the 2 Pz/?- 2 Pi /?. fine structu re line at 158 urn 
jHowe et al.ill99ll : lMookeriea et al.ll2003D . The radio obser- 
vations of carbon recombination lines (RRLs) provide a very 
good alternative way to trace the spatial distribution of ion- 
ized carbon. Classical C II regions which form part of high- 
excitation PDRs, are usually identified through the obser- 
vation of narrow (4 - 10 km s _1 ) ca rbon RRLs at frequen- 
cies > 1 GHz toward Hn regions iPankonin et al. 1 ll977l : 
Wvrowski et al.ll2000l) and have been well studied. The dif- 
fuse C II regions which form part of low-excitation PDRs, 
are identified through observations of carbon RRLs in ab- 
sorption at freque ncies below ~ 150 M Hz and in emission 
above ~ 200 MHz iPavne et al. Ill989l) . 



1 INTRODUCTION 

Molecular clouds are embedded in the diffuse interstel- 
lar medium (ISM) of the Galaxy where most of the 
gas is in an atomic form. The atomic gas is mostly lo- 
cated in two phases that are in pressure equilibrium, the 
cold neutral medium (CNM) with temperatures of T ~ 
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atomic gas is also associated with the denser molecu- 
lar c louds, possibly forming an extended halo or enve- 
lope iAndersson fc WannierlTl993l : iMoriartv-Schieven et alJ 



119971) . 

Far-ultraviolet (FUV) photons (6.0 < hv < 13.6 eV) 
from OB stars produce Photon Dominated Regions (PDRs) 
either at the interface between the H II region and the molec- 
ular cloud (classical high density PDRs) or in neutral com- 
ponents (atomic or molecular) of the diffuse interstellar 
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Diffuse C n regions were discovered more than two 
decades ago through the detection of carbon RRL in ab- 
sorption at 26.3 MHz (C631ct) in the directi on of Cas A 
llKonovalenko et al. I ll98Ct iBlake et al. Ill980h . Since then 
several RRLs spanning o ver 14 to 1 400 MH z hav e been ob- 
served in this direction (Ka ntharia et al~lll998l and refer- 
ences therein) and it remains the only diffuse C II region 
which has been studied and modelled in such detail using 
C RRL. A smooth transition of lines in absorption at fre- 
quencies below 115 MHz to lines in emission at frequen cies 
above 200 MHz has been detected JPavne et al. lll98St see 
Fig. 1). Moreover, this is the only direction in the Galaxy 
where high resolution (~ 1') image of t he diffuse C u re- 
gion in carbon RRLs near 332 MHz exist jKantharia et al~l 
Il998l . see Fig 2). Modeling the line emission in this direction 
shows that the carbon RRLs origi nate in cold regions wit h 
T e = 75 K and n e = 0.02 cm" 3 iKantharia et al. Ill99st) . 
Several molecular line studies have been performed towards 
the direction of Cas A. These observations do not detect 
much CO emission from the disk of Cas A, but identify 
mole cular clumps located primarily to the south east of the 
disk dTroland et alJll985l IWilson et aT1ll993l : iLiszt fc Lucas] 
1999). The diffuse Cn regions could be coexistent with the 
cold neutral m edium (CNM) that produce Hi absorption 
toward Cas A JPavne et al. Ill994l) or with the molecular 
comp onent of the interstellar medium (ISM) jErshov et al. I 
1987). In the direction of Cas A, a morphological comparison 
of C RRL distribution with that of Hi and 12 CO indicate 
that the C RRL emission is more likely associated with Hi 
gas iKantharia et al. "l ll99Sft . 

The upper state of the [C i] 3 Pi- 3 Po transition at 
492 GHz is only 23 K above ground. For this transition, 
the critical density for collisions with H2 molecules is only 
1000 cm" 3 JSchroder et al.lll99ll) . This implies that [Cl] is 
easily excited and also the line is easily detectable even when 
emitted by moderate density interstellar gas exposed only to 
a radiation field equal to the mean interstellar radiation field 
in the solar neighborhood. Thus [C 1] can be used as a re- 
liable tracer of the diffuse PDRs as well. Here we present 
[C 1] mapping observations at 492 GHz in the direction of 
Cas A in order to probe whether the [C 1] emission like the 
C RRL emission arise solely from the atomic CNM or the 
molecular phase also contributes to it. We have compared 
the [Cl] observations with the C270a, Hi, 12 CO 2-1 and 
13 CO 1-0 observations available in literature and have used 
PDR models to explain the neutral carbon and CO emission 
wherever the two spatially overlap. 



2 OBSERVATIONS AND DATA ANALYSIS 

We have mapped the Cas A region in the fine structure 
transition 3 Pi- 3 Po at 492 GHz of atomic carbon using the 
SubMillimeter Array Receiver for Two frequencies (SMART; 
Graf et al. 2002) on KOSMA (Winnewisser et al. 1986), a 
3-m submillimeter telescope located on the Gornergrat in 
Switzerland. SMART is a dual-frequency eight-pixel SIS- 
heterodyne receiver that observes simultaneously at 4 po- 
sitions (separated by 116") on the sky at two frequencies 
in the range 455-495 GHz and 795-882 GHz. The IF sig- 
nals are analyzed with array-acousto-optical spectrometers 
(array- AOSs). The array- AOS consists of 4 AOSs each with 
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Figure 1. [Cl] spectrum towards Cas A averaged over part of 
the observed region, limited in the west to Aa = —80" and in 
the north to A8 = 80". The Perseus arm and Orion arm compo- 
nents are clearly detectable. The effective velocity resolution of 
the spectrum is 0.63 km s" 1 . 



a bandwidth of 1 GHz and a spectral resolution of 1.5 MHz 
(Horn et al. 1999). This corresponds to a velocity resolution 
of 0.63 km s" 1 at the frequency of the [Cl] line. Owing to 
technical difficulties the higher frequency channel of SMART 
could not be used during our observations. Typical receiver 
noise temperature achieved at the center of the bandpass 
at 492 GHz is 150 K. Based on cross-scans on Jupiter at 
492 GHz we estimated the HPBW to be 55" and the beam 
efficiency (rj mb ) to be 50%. 

We observed a fully sampled map centred at Cas A 
(a 20 oo = 23 h 23 m 24 s ; <5 2 ooo = +58°48'.9), extending over 
~ 6'x 7'. The observations were done in the On-The-Fly 
(OTF) position switched mode, with the reference position 
being 2° south of the map centre. The map presented in 
this paper required 5 complete coverages of the region and 
the total integration time per position on the sky was 50 
seconds. The atmospheric calibration were done by measur- 
ing the atmospheric emission at the reference position to 
derive the opacity (Hiyama 1998) and the sideband imbal- 
ances were corrected for using standard atmospheric models 
(Cernicharo 1985). 

Data was analyzed using the GILDAS 1 spectroscopic 
data reduction package. Sinusoidal baselines were subtracted 
in order to correct for the standing waves. 



3 RESULTS 

3.1 Velocity Structure in [Cl] 

Figure shows the observed [C 1] spectrum averaged 
over part of the map restricted in the west to Aa = 
-80" and in the north to AS = 80". We identify 
mainly three strong emission features, at —47 km s" 1 and 
—38 km s" 1 corresponding to the Perseus arm and at 

1 http : //www . iram . f r/IRAMFR/GILDAS 
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Table 1. Results of Gaussian fits to the average spectral profiles 
of [Ci] 12 CO 2-1, 13 CO 1-0, C270a and optical depth spectrum 
of Hi. 



Line 


Component 


Ccntroid 
km s _1 


Width 
km s _1 


[Ci] 


Perseus arm 
Orion arm 


-47.4 ±0.1 
-1.7 ±0.1 


2.9 ±0.2 
3.1 ±0.4 


12 CO 2-1 


Perseus arm 
Orion arm 


-46.9 ±0.01 
-40.3 ± 0.01 
-36.2 ± 0.03 
-1.1 ±0.01 


3.0 ±0.05 
4.5 ±0.07 

3.0 ±0.07 

2.1 ±0.02 


13 CO 1-0 


Perseus arm 
Orion arm 


-46.8 ±0.05 
-39.7 ±0.15 
-1.1 ±0.06 


2.1 ±0.1 
6.1 ±0.3 
1.8 ±0.14 


C270a 


Perseus arm 


-46.8 ± 0.6 


5.1 ±0.9 
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Figure 2. Left panel: Average spectra of [Ci], 12 CO 2-1, 13 CO 
1-0 emission. Right panel: Bottom: Average spectrum of C270o 
in units of line to continuum ratio and Top: Average spectrum of 
H I optical depth. All spectra are averaged over the disk of Cas A. 



— 1 km s -1 corresponding to the local Orion arm. Since the 
—38 km s _1 feature and the Orion arm feature do not ap- 
pear at all positions within the mapped region, the spectrum 
obtained by averaging over the entire mapped region does 
not show them clearly. Thus, in order to provide a fair rep- 
resentation of all velocity components detected in the region 
we show here the average spectrum over the selected region. 

Figure |21 shows a comparison of the emission spectra of 
[Ci], 13 CO 1-0, 12 CO 2-1, C270a (at 332 MHz) and the 
optical depth spectrum of Hi, all averaged over the disk 
of Cas A. The CO spectra are from the observations by 
[ Liszt fc Lucasl Jl99Sl>. t he C270a spectrum was observed by 
Kantharia et al. (1998), and the H I optical depth spectrum 
is from lSchwarz et al] 1119971) . Tabled presents the results of 
fitting Gaussian components to the average emission spectra 
of the different tracers. The —47 km s _1 Perseus arm and the 

— 1 km s _1 Orion arm features (not in C270a) are detected 
most clearly in all the emission and absorption spectra. The 




Figure 3. Distribution of [Ci] emission from Cas A. The inten- 
sities are integrated between —50.3 and —44.1 km s _1 . Contours 
range 28 to 98% (in steps of 10%) of the peak (9.7 K km s" 1 ). 
Circle denotes the extent of the Cas A disk. 



CO 2-1 emission shows two more Perseus arm features 
nominally around -39 and -36 kms -1 . The [Ci], 13 CO 
and Hi data do not resolve the two Perseus arm features 
at —39 and —36 km s _1 , rather show a single feature at 
~ —38 km s _1 . Although the C270a spectrum shows mainly 
the —47 km s _1 feature, the same gas as traced by C RRLs 
with different as also shows the —38 km s" 1 feature bo th in 
emission and absorption (c/. Fig. l. lPavne et al. lll98SlL 

All tracers appear to show rather similar velocity com- 
ponents. Thus based on velocity features it is not possible to 
isolate the contribution of atomic and molecular components 
of the neutral ISM towards the observed [C i] emission. 

Rest of the paper discusses the —47 km s _1 Perseus arm 
feature which is detected in all tracers. 



3.2 [C i] emission from the Perseus Arm at 

-47 km s" 1 

Figure [3] shows the intensity map of [C i] emission from the 
Perseus arm clouds along the line of sight to Cas A (the 
circle denotes the extent of Cas A), integrated between the 
velocities of —50.3 and —44.1 km s _1 . The observed [Ci] 
emission toward Cas A consists of two distinct morphologi- 
cal features. The first feature is the bright emission detected 
towards the south-east and the second is the filamentary 
structure to the west of the map extending almost north- 
south. There is almost no emission detected from the north- 
ern part of the disk of Cas A. 

Figure 2] shows the channel maps of [C i] emission be- 
tween velocities —49.7 km s _1 to —44.7 km s _1 . It shows 
that the emission from the cloud to the south-east covers the 
largest range of velocities between —48.5 to —45.9 km s" 1 . 



4 B. Mookerjea, N. G. Kantharia, D. Anish Roshi & M. Masur 



0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 0.5 1.0 1.5 2.0 



58 53 
52 
51 
50 
49 
48 
47 
46 
45 
58 53 
52 
51 
50 
49 
48 
d 47 

LU 

Q 46 
45 
58 53 
52 
51 
50 
49 
48 
47 
46 
45 



O 
i= 
< 



-49^1 KM/S 




-48.5 KM/S 













i 


<4> o ( 






liS-mliir 




Figure 5. Comparison of [Cl] (gray scale) channel (of 
1.4 km s width) map w ith C270» optical depth (black contours) 
jKantharia et al71ll998l) centred at a veloc ity of —48.0 km s —1 . 
The circles show the CO clouds identified bv lWilson et al.Ul993t) . 
The contour levels are 1, 2, 3, 4, 5, 7, 9 and 11 in units of 0.001. 
The greyscale ranges between 2 and 9.8 K km s — 1 . 
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Figure 4. Channel maps of CI emission. The circle denotes the 
extent of Cas A. Velocity is labelled in each panel. The grey scale 
ranges from 0.32 to 2 K and contours are plotted for 0.48, 0.68, 
0.96 and 1.35 K. Note the wider south-east component of [Cl]. 



The linear filamentary extension in the west in contrast is 
detected only between —49.1 to —47.2 km s _1 . 

The morphology of [C i] emission rules out interaction 
of this gas with the SNR. As shown in Fig. 2] the [C i] emis- 
sion is seen mostly from the southern part of the supernova 
remnant; the easte rn part coincides wi th the molecular cloud 
seen in that region dWilson et alll993l) which is also believed 
to lie in the Perseus arm but not interacting with Cas A. If 
there was interaction between [C i] line emitting region and 
Cas A then it would lead to one of the following being true: 
a) if the line emission is affected by the continuum from 
Cas A then it should show good morphological correlation 
with the supernova remnant (SNR). However this is not the 
case here, b) peculiar spectral features from the regions close 
to Cas A e.g. broad lines or stronger lines. No such trend 
is seen and this also helps rule out any interaction between 
the two. Thus, we believe that the [C i] gas is not interacting 
with the supernova remnant. 



3.3 Comparison of [Cl], C270a & Hi distribution: 
Atomic Gas 

Figure|^|shows the [C i] emission superposed on the contours 
of C270a optical depth jKantharia et al"lll998l) both at ve- 
locities of —48 km s . While the [C i] emission overlaps par- 
tially with the C270a emission, it also appears to surround 



the C270a emitting region and there are some significant dif- 
ferences in morphology between the two. The C270cv. peak 
optical depth is observed from the central part of Cas A and 
is not coincident with the [C i] peak. The [C i] peak is located 
to the south of the C270a peak which possibly indicates the 
chemical stratification from C + through C° to CO. 

The location of the peak of the C RRL emission at 
the centre of the SNR in Cas A is indeed intriguing. As- 
suming that the observed line width of recombination lines 
at 34.5 MHz was ent irely due to radiation broadening, 
jKantharia et al.1 lll998l) put a lower limit on the distance of 
115 pc between the line forming region and Cas A indirectly 
ruling out t he association of thi s gas with the SNR. On the 
other hand. lKassim et alJ lll995l) have reported flatter spec- 
trum in the central arcmin or so of Cas A at frequencies 
below 330 MHz. The common location of the recombination 
line gas and the absorbing g as is intriguing and argues for 
a common origin. However, iKassim et al.l lll995l) argue that 
the thermal absorbing gas is located within Cas A in which 
case a common origin of the gas is ruled out. Moreover the 
carbon recombination line forming gas has a continuum op- 
tical depth of 6.4xl0 -4 at 74 MHz wh ich is much less than 
the optical depth (1.3) of the gas that lKassim et all dl995l) 
detect. We conclude that the positional coincidence of the 
carbon recombination line forming gas and the flat spec- 
trum region towards Cas A does looks fairly intriguing and 
demands further investigation. 

The line of sight to Cas A has provided a 
uniquely detailed set of observations of RRLs from 
highly excited states of singly ionized carbon (e.g. 
Konovalenko et al. ll98Ct lErshov" t al. Il987t iPavne et al. I 



19891: lAnantharamaiah et al. Hl994HKantharia et al. I l998t. 
Interpretation of these RRLs using models considering 
dielectronic-like recombination suggest that with the excep- 
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tion of the recombination line width, all of the Cas A carbon 
recombination line and A21 cm Hi absorption line data can 
be attributed to a region where the physical c onditions are 
^typica l of the cold neutral medium of the ISM jPavne et al. I 

We also compar ed the [C il emission with the Hi distri- 
bution across Cas A dSchwarz et al Jl997l) . H 1 is observed in 
absorption across the entire Cas A disk and the morphology 
varies with the radial velocity with maximum optical depth 
being observed around —48.2 km s" 1 . The peak absorption 
at this velocity arises in the south-east region which roughly 
coincides with the peaks seen in [C i] at similar velocity. H I 
absorption is also seen from the western half of Cas A with 
the peak shifted slightly towards the south. A similar peak 
is also observed in [C i] emission. 

The [C i] distribution is thus partly similar to both H I 
and C270a morphology. We compare the [C i] emission with 
the molecular component of the ISM traced by the CO emis- 
sion in the next subsection. 

3.4 Comparison of [C i] and molecular line 
emission toward Cas A 

The left and middle panels of Fig. |S| show the distribu- 
tion of 13 CO 1-0 and 12 CO 2-1 intensities (in contours) 
integrated between v elocities of —50.3 and —44.1 km s~ 
iLiszt fc Lucaslll999l) . superposed with the [C i] distribu- 
tion (in grey). While the original resolution of the 13 CO 
1-0 map is 60", the 12 C O 2-1 data h a s been smoothed to 
60". In an earlier paper IWilson et al] lll993f) observed the 
CO (and 13 CO) emission from the Perseus arm clouds in 
the direction of Cas A, with an ang ular re s olutio n ~ 22". 
The CO clouds identified by IWilson et al] lll993l) are also 
marked by circ l es in Fig. |S] We note that the CO data from 
iLiszt fc Lucas ] <|l999T) do not sh ow the clumped structure 
observed bv I Wilson et all lll993|) and we attribute it to the 
difference in resolution between the two datasets. 

Figure |S| shows that the [C i] emission is detected from 
most of the CO emitting clouds and the [C i] peaks lie at the 
outer boundaries of the clumps. Substantial [C i] emission is 
detected with no corresponding CO emission from the dif- 
fuse PDRs to the north. The extreme right panel of Fig. HJ 
shows the integrated [ C il intensities over plotted with the 
CO clouds identified bv lWilson et alJ Il993t) . It clearly shows 
that although [C i] emission is detected from most of the 
clouds, the peaks of [C i] emission tend to lie at the bound- 
aries of the CO clouds. We point out here that the CO clouds 
detected by I Wilson et al] il993T) towards the western edge 
of the map belong to the —36 km s" 1 Perseus arm feature, 
in contrast to the —47 km s" 1 feature in [Ci] emission be- 
ing discussed in this paper. The relative position of the [C i] 
and CO emission peaks is consistent with the clumpy PDR 
scenario, in which PDRs form on the surface of the molecu- 
lar clumps that are embedded in an interclump medium of 
lower density in which CO is dissociated to atomic carbon. 

From the above comparisons it is clear that [C i] traces 
both CO and cold Hi morphologies partially and is consis- 
tent with the PDR scenario when compared with the C + re- 
gions traced by the carbon recombination lines (e.g. C270q) 
and the molecular clouds traced by CO. We note here, that 
absorption studies of C RRLs necessarily imply that the 
emitting gas lies in front of the strong continuum emitting 



SNR, no such restriction is however applicable to either the 
CO or the [Ci] lines. Thus, beyond the similarities in the 
emission velocities of [C i] CO and the C270a and the mor- 
phological similarity to a PDR-like structure, the present ob- 
servations cannot conclusively rule out the possibility that 
the molecular (and neutral) and the ionized gas may not be 
occupyng the same volume of gas. 



4 NEUTRAL CARBON AND CO COLUMN 
DENSITIES 

We have estimated the C° and CO column densities based 
on LTE approximation. For these calculations we have used 
the observed [C i] and 13 CO 1-0 line intensities and an exci- 
tation temperature, T ox , of 20 K. Further we have assumed 
both [Ci] and 13 CO 1-0 to be optically thin. The assumed 
T cx of 20 K is justified because based on the 12 CO and 
13 CO observations the kinetic temperatures of the molecu- 
lar clumps in the Perseus arm along the direction of Cas A 
is ~ 20 K JWilson et al.lll993h . 

In the Perseus arm clouds with velocity ~ —47 km s" 1 
we estimate the column density of neutral carbon, N(C), to 
be between 2xl0 16 cm" 2 and 1.3xl0 17 cm" 2 . Near the [Ci] 
peaks in the south-east and also to the west, N(C) is larger 
than 7 10 16 cm" 2 , while to the north and in the periphery 
of the [C i] emitting shell the column densities are close to 
3 10 16 cm" 2 . 

Within the region mapped by us in [Ci], the esti- 
mated 13 CO column density (in a beam of 1') varies be- 
tween 5xl0 14 cm" 2 and 9x 10 15 cm" 2 . This corresponds 
to N(CO) between 3.0 x 10 16 cm" 2 and 5x 10 17 cm" 2 , as- 
suming [ 12 CO/ 13 CO]=60 . The column density o f 13 CO de- 
rived from the 1' dataset of ILiszt fc^ucas] (^99l ) is sli ghtly 
lower than the values obtained bv Wi lson et al] (1993), but 
matches reason ably well with the p revious l.'l 13 CO 1-0 
observations bv lTroland et all <ll985l) . 

Over the limited region in which the 13 CO and [C i] 
emissions overlap, the C/CO abundance ratio varies between 
0.08 and 1.8. At the positions of the CO clouds F and G the 
C/CO abundance ratio is ~ 0. 2, a value typically found at 
the c entre of dense cloud cores llGerin et aljll998l : IOka et al] 
2004). In the periphery of the CO clouds the value rises to 
1.7 at positions where the [Ci] emission peaks : this value is 
typical of the diffus e, translucent clouds lllngalls et al ll997t 
iBensch et al"ll2003ft . 



5 PDR INTERPRETATION OF THE 
OBSERVED [Ci] EMISSION 

The observed morphology of [C i] though not entirely cor- 
related with the emission from the atomic or the molecular 
phases of the ISM, is consistent with the phenomenological 
understanding of PDR structure. [C i] is typically ubiquitous 
and towards the direction of Cas A, it appears to be aris- 
ing from both the molecular clumps as well as the diffuse 
PDR. Most importantly, the [C i] emitting regions show no 
evidence for interaction with the SNR in Cas A. Here we 
explain the observed [C i] emission in terms of theoretical 
m odels for PDRs. We ha ve used the plane-parallel models 
bv lKaufman et al. I (Il999l) for this purpose. The PDR model 
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Figure 6. Integrat ed intensity distribu tion of [Cl] (Grey) overlayed with integrated intensities of 13 CO 1-0 (left panel) and 12 CO 
2-1 (middle panel) jLiszt & Lucadll999l) . The 13 CO 1-0 data has a resolution of 60". The 12 CO 2-1 dataset has been smoothed to a 
resolution of 60" for comparison. The large circle denotes the angular extent of Cas A. Extreme right panel shows both in grayscale 
and contours the in tegrated intensity of [Cl]. The smaller circles in all panels show the position and sizes of CO clouds identified by 
IWilson et"aT]lll993ft . 



used here takes into account a detailed chemical network, 
radiative transfer and thermal balance and considers a 1- 
dimensional slab of material exposed to a FUV radiation 
field in a face-on geometry. The plane-parallel PDR model 
calculates the intensity of various transitions of different 
chemical species as functions of two parameters, the local hy- 
droge n density ( n ppR,) and the incident FUV radiation field. 
iKaufman et al. I <ll999l) have created a database of PDR 
models, for local hydrogen densities hpdr, between 10 & 
10 r cm" 3 and FUV radiation field between 0.3 to 3x 10 6 Go, 
where Go is the FUV flux measured in units of the "Habing 
field", which is taken to be 1.6xl0 -3 erg cm -2 s _1 sr _1 
<Habindll96sl) . 

In the absence of any nearby stars the far-ultraviolet 
(FUV) radiation field illuminating the [C i] emitting neutral 
atomic and/or molecular phase of the ISM will be of the 
order of 1-5 Go. For the following PDR analysis we assume 
the FUV flux to be ~ 1 G . 



5.1 Hydrogen densities 

We have used the observed CI/ 12 C02-1 intensity ratios to 
derive an independent estimate of the hydrogen volume den- 
sities in the [C i] emitting regions using the PDR models. For 
convenience from now on we split our discussion into two 
parts, the first deals with the south-eastern [C i] emitting 
regions which are largely concomitant with the CO emission 
and the second deals with the north-south stretching fila- 
mentary emission to the west of the map. The aim of this 
separation is also to look for possible differences between the 
two parts. 

Figure □ shows the distribution of the CI/ 12 C02-1 in- 
tensity ratios in regions where both intensities are more than 
three times the respective noise levels. Most of the pixels in 
both the south-eastern cloud and the western filament show 
intensity ratios between 0.5 and 1.0. At the position of the 
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Figure 7. [Cl]/ 12 C02-l intensity ratio towards Cas A. Only 
positions with both [Cl] and 12 CO 2—1 intensities larger than the 
3ct values for the respective maps are included. The contours of 
[Cl] emission are overla yed for comparison. Also marked are the 
CO clouds identified bv lWilson et alj 1 1993). The map is centred 
at the nominal centre being used throughout the paper, viz., 02000 
= 23 ,l 23 m 24 3 ; feooo = +58°48'.9. 



[Cl] emission peak to the south-east the ratio is 1.3 and 
towards the northern edge of the south-eastern cloud the 
ratio rises from 1.8 to 4. In the western filament also the 
value rises to as high as 4 towards north. Clouds F, G and 
H show ratios of ~ 0.5, while cloud E has a value of ~ 1.3. 

Figure |H| shows the variation of the CI/ 12 C02-1 inten- 
sity ratios as a function of the UV radiation field and the 
loc al hydrogen den s ities a s estimated by the PDR models 
bv lKaufman et al. I Jl99fih . For convenience we have drawn 
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Figure 8. PDR model calculations of variation of CI/ 12 C02-1 
intensity ratios as a function of hydrogen densities and incident 
UV radiation field iKaufman et al. Il999t) . The contours (marked 
with their values) depict some of the typical ratios found toward 
Cas A. The horizontal line corresponds to the UV radiation field 
of 1 G . 



contours corresponding to some of the typical ratios found in 
the direction of Cas A. We note that according to the PDR 
models, the observed CI/ 12 C02-1 intensity ratio constrains 
the local hydrogen density almost independent of the UV 
radiation field assumed. CI/ 12 C02-1 intensity ratio of 4 or 
more corresponds to densities less than 100 cm -3 . The most 
commonly seen ratio of 0.5 suggests densities of the order of 
10 3 ' 5 cm -3 . Thus the hydrogen density of the [Ci] and CO 
emitting gas varies from 100 to 3000 cm -3 . 

Since the [C i] emission only partially overlaps with the 
carbon RRL emission as well as the CO emission in a man- 
ner apparently consistent with the PDR scenario, it is in- 
teresting to compare the neutral hydrogen densities in the 
different regions. Based on latest models of carbon RRL 
emission, the atomic hydrogen densities are estimated to be 
around 150 cm~ 3 for a n electron temperature, T e , of 75 K 
llKantharia et al"lll99ct) . CO observations suggest hydrogen 
volum e densities of ~ 10 3 cm -3 with localized higher den- 
sities JWilson et al.lll993D . We find that the southern parts 
of the mapped region [C i] emission arises primarily from re- 
gions in which hydrogen is mostly molecular with average 
densities of 10 3 cm' 3 . To the north, the densities are lower 
and the [C i] emitting region has hydrogen primarily in the 
atomic phase. This is shown by the higher [C l]/ 12 C02-l in- 
tensity ratios and the lower densities derived from the car- 
bon RRL emission which peaks in this region and is known 
to arise in the atomic medium. The [C i] emission thus not 
only morphologically but also in terms of the continuity of 
density parameters acts as the transition phase between the 
atomic and the molecular phase of the neutral ISM. 



5.2 Volume and Area Filling factors of the [C i] 
emission 

In order to quantify our results better we have selected a few 
representative positions and derived densities at those posi- 
tions for an UV field of x — 1 Go. This analysis also aims 
to characterize the clumpy nature of the emitting regions 
along the lines of sight towards Cas A. The positions cho- 
sen correspond to t he centres of the CO clouds E, F and G 
JWilson et al.ll 993) from which [C i] emission was clearly de- 
tected. In addition, we have chosen the primary [C i] peak in 
the entire mapped region, lying in the south-eastern molec- 
ular cloud and the secondary peak found to the north of the 
filamentary emission to the west. Finally we have chosen a 
representative position with average emission characteristics 
in the western filament. 

Table^summarizes details of the selected positions and 
several parameters observed and derived from PDR models 
and obtained from literature. 

As discussed in Sect. 15. li the local hydrogen volume den- 
sities, npDR, at the selected positions are estimated from the 
observed CI/ 12 C02-1 intensity ratios using the PDR mod- 
els (Col. 7, Table 0. We find n PDR for the CO clouds F 
and G to be ~ 10 3 cm -3 , while cloud E has a much lower 
density of only 300 cm -3 . These densities are in agreement 
with the values derived bv I Wilson et al.l <ll993h using much 
higher resolution (22") observations (Col. 12 Table 0. At 
the position of the two [Ci] peaks to the south-east and in 
the western filament the densities are ~ 500 cm -3 , while for 
most of the positions on the western filament the density is 
~ 10 3 cm -3 . This implies that the high [Ci] intensities arise 
from the more diffuse interclump medium. 

The N(H 2 ) per 60" beam (Col. 5, Table for the 
individual positions have been estimated as described in 
Sec. H For the clouds E & G the N(H 2 ) per beam (60") 
are lower by more than a factor of 3 than the N(U" 2 ) mea - 
sured per 22" (Col. 11, Tabled bv IWilson et alJ Il993l) . 
while for cloud F the two numbers agree. At those posi- 
tions where the 13 CO 1-0 emission lies below the noise lim- 
its of the dataset, we have assumed the noise level to derive 
N(H2)< 5.6 x 10 20 cm -2 . Using the hydrogen column densi- 
ties and assuming that the cloud extends the same distance 
along the line of sight as the linear size of the beam we have 
derived the beam-averaged volume densities (n av ; Col. 6, 
Table HJ. We note here that n(H 2 ) w ii SO n (Col. 12, Table |5J 
has also been derived from the N(H2) w iison in the same way. 

The volume filling factors of the emitting clouds can 
be estimated as <j>v=npr>n/n a . v . For most of the positions 
except for Cloud G and the average position in the west- 
ern filament the volume filling factors are higher than 50% 
(Col. 10, Table This is in contrast to the typical volume 
filling f actors of 5% as foun d in the Galactic star forming 
regions jKramer et aljl2004h . The observed high volume fill- 
ing factors thus indicate the rather diffuse nature of the [C i] 
emitting regions. 

Comparison of the absolute [C i] intensities, Ici,th (Col- 
umn 7, Table 0, estimated by the PDR models which best 
reproduce the observed CI/ 12 C02-1 ratios, with the ob- 
served [C i] intensities (Ici,ob s ; Col. 4, Table |5J provides in- 
formation about the area filling factors (Col. 9, Tabled of 
the [C i] emitting clouds. The area filling factor is defined as: 
</>A=Ici,obs/Ici,th. At the selected positions in the direction 
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Table 2. Parameters observed and derived from PDR models at selected positions in the Perseus arm. nprjR is local hydrogen volume 
density derived by comparing the observed CI 12 C02-1 intensity ratios with predictions of the PDR models. N(H2) is the H2 column 
density per 60" beam, cfi\ is the area filling factor of CI e mission defined as Ig ba/Ici th an d 4>V^ S the volume filling factor defined as 
"pdh/%v Quantities with the subscripts wilson are from Wils on et alj §1993). 



(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 

Name Act AS I C i,obs N(H 2 ) n av n PDR I CIjth </> A </> v N(H 2 ) wi i son n(H 2 ) w ji, 

" " K km s" 1 10 20 cm -2 cm -3 cm" 3 K km s" 1 10 20 cm -2 cm -3 



Cloud E 58. 0. 5.0 5.6 180 300 10.9 0.5 0.6 20 <800 

Cloud F 58. -87. 5.3 39 1200 2000 14.1 0.4 0.6 40 1000 

Cloud G 87. -174. 6.7 52 1600 5000 13.0 0.5 0.3 90 2000 



CI Peak West -87. 29. 2.9 5.6 180 400 10.9 0.3 0.5 

CI Peak SE 29. -87. 8.1 18 560 630 12.6 0.6 0.9 



CI West Avg -87. -58. 4.8 11 340 1778 14.1 0.3 0.2 



towards Cas A the [Cl] area filling factor is on an average 
equal to 0.5. 



6 SUMMARY 

We have used the [Cl] emission at 492 GHz to probe the 
PDRs seen in the direction of the supernova remnant Cas A. 
We detect [Cl] emission at velocities around —47 and at 
—39 km s _1 from the Perseus arm clouds, and at —1 km s _1 
from the local Orion arm cloud. The [C 1] emission does not 
show any strong morphological correlation with the contin- 
uum emission from Cas A. Further, we do not detect any 
obvious broadening or strengthening of spectral lines closer 
to the continuum peak. Both of these conclusively rule out 
that the [C 1] emitting gas is interacting with the SNR in 
Cas A. 

We detect [C 1] emission from both the diffuse CNM (as 
traced by the C RRLs) and the denser molecular cloud. Mor- 
phological comparison of the [C 1] emission with the carbon 
RRL (C270q), Hi optical depth and CO distribution appear 
to be consistent with a PDR scenario in which the dominant 
carbon-bearing species changes from C + in the diffuse region 
through C° to CO in the more dense molecular clumps. We 
note with caution that despite the morphological and kine- 
matic evidences the present observations cannot rule out the 
possibility that the [C 1] and C RRL emitting regions may 
not completely occupy the same volume of space. 

We estimate the C° column density, N(C), to be be- 
tween 2 10 16 and 1.3 10 17 cm -2 . At the position of the pre- 
viously identified CO clouds we estimate the C/CO abun- 
dance ratio to be around 0.2, similar to the typical value 
found in dense cloud cores. Most of the [C 1] emitting re- 
gion however shows C/CO abundance ratios > 1, typical of 
diffuse, translucent clouds. 

We have derived hydrogen volume densities of the [C 1] 
emitting regions from the [Cl]/ 12 C02-l intensity ratios us- 
ing PDR models. We find that the local hydrogen densities 
vary from ~ 10 2 cm -3 near the [Cl] peaks to 10 3 ' 5 cm -3 
around CO peaks. This also shows a continuity in the den- 
sity structure in the region between the the diffuse C 11 re- 
gion and the denser molecular clouds. We estimate the area 
filling factor of the [C 1] emitting regions to be ~ 0.3, while 
the volume filling factor is on an average > 30%. This fur- 



ther suggests that most of the [C 1] emission stems from the 
diffuse PDR also traced by the carbon recombination lines. 



ACKNOWLEDGMENTS 

We thank the referee, Jay Lockman, for comments which im- 
proved the clarity and presentation of the paper. We thank 
H. Liszt for allowing us to use the CO datasets. This material 
is based upon work supported by the Deutsche Forschungs 
Gemeinschaft (DFG) via grant SFB494 and the National 
Science Foundation under Grant No. AST-0228974. This re- 
search has made use of NASA's Astrophysics Data System. 



REFERENCES 

Anantharamaiah, K. R., Erickson, W. C, Payne, H. E., 
Kantharia, N. C, 1994, ApJ, 430, 682. 

Andersson, B.-G., & Wannier, P. G. 1993, ApJ, 402, 585 

Bensch, F., Leuenhagen, U., Stutzki, J., & Schieder, R. 
2003, ApJ, 591, 1013 

Blake, D. H., Crutcher, R. M., Watson, W. D., 1980, Na- 
ture, 287, 707 

Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 
Ershov, A. A., Lekht, E. E., Smirnov, G. T., Sorochenko, 

R. L. 1987, Sov. Astron. Lett., 13, 8 
Gerin, M., Phillips, T. G., Keene, J., Betz, A. L., & Boreiko, 

R. T. 1998, ApJ, 500, 329 
Habing, H. J. 1968, Bull. Astr. Inst. Netherlands, 19, 421 
Hollenbach, D. J., Tielens, A. G. G. M. 1999, Rev of Mod 

Phys., 71, 173. 

Howe, J. E., Jaffe, D. T., Genzel, R., & Stacey, G. J. 1991, 
ApJ, 373, 158 

Ingalls, J. G., Chamberlin, R. A., Bania, T. M., et al, 1997, 
ApJ, 479, 296 

Kantharia, N. G., Anantharamaiah, K. R., Payne, H. E., 
1998, ApJ, 506, 758 

Kaufman, M. J., Wolfire, M. G., Hollenbach, D. J., Luh- 
man, M. L., 1999, ApJ, 527, 795 

Kassim, N. E., Perley, R. A., Dwarakanath, K. S., & Er- 
ickson, W. C. 1995, ApJ, 455, L59 

Konovalenko, A. A., Sodin, L. G., 1980, Nature, 283, 360 

Kramer, C, Jakob, H., Mookerjea, B., et al., 2004, A&A, 
424, 887 



[Cl] observations toward Cas A 9 

Kulkarni, S. R., & Heiles, C. 1987, ASSL Vol. 134: Inter- 
stellar Processes, 87 

Liszt, H., & Lucas, R. 1999, A&A, 347, 258 

McKee, C. F., Hollenbach, D., & Wolfire, M. G. 2004, The 
Dense Interstellar Medium in Galaxies, 395 

Mookerjea, B., Ghosh, S. K., Kaneda, et al., 2003, A&A, 
404, 569 

Moriarty-Schieven, G. H., Andersson, B.-G., & Wannier, 

R G. 1997, ApJ, 475, 642 
Oka, T., Iwata, M., Maezawa, H., et al., 2004, ApJ, 602, 80 
Pankonin, V., Barsuhn, J., Thomasson, P., 1977, A&A, 54, 

335 

Payne, H. E., Anantharamaiah, K. R., Erickson, W. C, 

1989, ApJ, 341, 890 
Payne, H. E., Anantharamaiah, K. R., Erickson, W. C, 

1994, ApJ, 430, 690 
Schwarz, U. J., Goss, W. M., & Kalberla, P. M. W. 1997, 

A&AS, 123, 43 
Schroder, K., Staemmler, V., Smith, M. D., Flower, D. R., 

& Jaquet, R. 1991, Journal of Physics B Atomic Molecular 

Physics, 24, 2487 
Troland, T. H., Crutcher, R. M., & Heiles, C. 1985, ApJ, 

298, 808 

Wilson, T. L., Mauersberger, R., Muders, D., Przewodnik, 

A., & Olano, C. A. 1993, A&A, 280, 221 
Wolfire, M. G., McKee, C. F., Hollenbach, D., & Tielens, 

A. G. G. M. 2003, ApJ, 587, 278 
Wyrowski, F., Walmsley, C. M., Goss, W. M., & Tielens, 

A. G. G. M. 2000, ApJ, 543, 245 



